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Abstract. We investigate a possible explanation for the NuTeV anomaly in terms of a nuclear 
correction difference between m,, and dy distributions. Analyzing deep elastic scattering and Drell- 
Yan data for nuclear targets, we try to determine the correction difference and its effect on the 
anomaly. We find that the difference cannot be precisely determined at this stage due to the lack 
of data which are sensitive to the difference. Therefore, it is difficult to draw a solid conclusion 
about its effect although the anomaly could be explained, at least partially, by this kind of nuclear 
correction. 



INTRODUCTION 

Weak-mixing angle is one of the fundamental constants in the standard model. Us- 
ing neutrino and antineutrino deep inelastic scattering data, the NuTeV collaboration 
found that their angle ^ir? dw = 0.2277 ± 0.0013 (stat) ± 0.0009 (syst) |1] in the on- 
shell scheme is significantly larger than the average of other measurements, sin 6w = 
0.2227 ±0.0004. This difference is called "NuTeV anomaly". Before discussing any ex- 
otic explanations, we need to investigate possible sources from nucleonic and nuclear 
structure. For example, the strange asymmetry s — s and QED effect contribute to the 
difference. In addition, appropriate nuclear corrections should be taken into account be- 
cause the NuTeV target is the iron. Among the nuclear corrections, it was pointed out 
that the nuclear modification difference between and dy distributions could affect the 
sin Ow determination [2]. Then, this nuclear modification is determined by using world 
nuclear-correction data on the structure function F2 and the Drell-Yan cross section in 
order to find the possible origin of the NuTeV anomaly 0]. This analysis was done 
by using a analysis technique developed for determining nuclear parton distribution 
functions L^. The following discussions are based mainly on the works in Ref. 



GLOBAL ANALYSIS FOR DETERMINING NUCLEAR 
CORRECTION DIFFERENCE BETWEEN Uy AND Dy 

From the neutrino and antineutrino deep inelastic scattering data, it is possible to obtain 
the weak mixing angle by the Paschos-Wolfenstein relation, R = (c^^ — Ojvc ) / (*^cc ~ 



^cc) ~ V2 — sin^ Ow, where a^^ and a^^ are charged-current (CC) and neutral-current 
(NC) cross sections, respectively. Various correction factors to this relation are discussed 
in Ref. [2] . Among them, we investigate a possible nuclear correction difference between 
Uv and dy. These valence-quark distributions should satisfy the baryon-number and 
charge conservations so that they have certain restrictions [2]. In principle, they could 
be determined by a global analysis of world data on nuclear structure functions. 
Valence-quark distributions are defined in nuclei as |0] 

uf;{x,Q }=WuAx,Q ,A,Z) J , 

"v [x,Q)= Wd,[^, Q ,A,Z) , (1) 

where Wu^, and w^^, indicate nuclear modification factors. It should be noted that the 
modification factors w/ are used at any points in this work |2. 0] although they are 
defined only at Q^=\ GeV^ in Ref. [4]. The difference is defined Awy = W;,,, — w^,,, and 
it is expressed by four parameters, a'„ h[„ c[„ and d[, at Q^=\ GeV^ (= Q^): 

where the value in Ref. |4] is used for jSy. The parameters are determined by a global 
analysis by minimizing the total x^: = Y,j{Rj"'" — R'-"'")^ / {Oj"'")^ , where Rj indi- 
cates the structure-function and Drell-Yan cross-section ratios, F^/F^' and of^/o^ . 
The experimental error (aj'"'")^ is given by the quadratic summation of statistical and 
systematic errors. Uncertainty of the obtained Awy is calculated by the Hessian method. 



RESULTS 

The nuclear valence-quark distributions are de- 
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fined at Q =1 GeV with the parameters in Eq. Q'-icev' 
([T]). We use the antiquark and gluon distributions 
which are determined in |4]. The nucleonic dis- | 
tributions on the right-hand sides of Eq. ([T]) are 
taken from the MRSTOl distributions. The distri- S .ooi 
butions are evolved to the experimental points of 

F^/Ff and o^y /^dy by the DGLAP evolution ^ 

equations. Then, the parameters are determined °°' " ' 

in comparison with the data so as to minimize ^^r^TTT,.. , ^, • j a 
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1, and the uncertainty of Aw,, is shown by the band. We find that the difference Awy 
is a small quantity. In fact, the nuclear modifications for and d^ are assumed to be 
the same in most theoretical calculations of nuclear structure functions. The difference 



could be related to QED effects on the nuclear parton distributions. However, its precise 
determination is not possible from experimental data because the uncertainty is an order 
of magnitude larger than the obtained distribution in Fig. 1. 

Because the Paschos-Wolfenstein relation was not directly used in the NuTeV analy- 
sis, we need to take a weighted average over the NuTeV kinematics for discussing the 
anomaly. For example, there are few data in the large-x region. Fortunately, such a calcu- 
lation method is provided in Ref . [ 1 ] . Because the NuTeV definition of the nuclear parton 
distributions is slightly different from ours, their relations should be found. According to 
the NuTeV convention, the nuclear distributions are defined as xu^ = {Zu*p+ Nu*i^) /A 

and xd^ = {Zd*p +Nd*^) /A, where the NuTeV distributions are denoted with the asterisk 
(*). Comparing these expressions with Eq. (HI), we find that the modification difference 
Awy corresponds to isospin-violating distributions of the NuTeV collaboration: 

5 M* = ulp - d*„ = +Awy xuy, 5d* = d*p - u*,^ = -Aw^ x dy . (3) 

Information is given in Ref. [1] for calculating these isospin-violating distributions on 
the weak-mixing angle: 

A(sin2 0iy) = - J dx{F[5u*„x] 5u*,{x) + F[5d;,x] 5d*{x) }, (4) 

by using the provided functionals F[5ul,x] and F[5d*,x]. 

Using the distribution Awv together with Eqs. © and ©, we obtain A(sin^ Bw) — 

0. 0004 ±0.0015. It should be noted that the distribution Aw,, is calculated at 2^=20 
GeV^, which is approximately the average value of the NuTeV experiment. The 
magnitude is not large enough to explain the anomaly (0.0050). However, it does not 
mean that the considered nuclear modification mechanism should be ruled out. There 
is a possibility that the error is underestimated because the error correlation effects 
with antiquark and gluon distributions are not taken into account. In any case, it is 
obvious from Fig. 1 that Awy cannot be determined at this stage, so that a precise 
numerical estimation is not possible. Future efforts are needed for clarifying the nuclear 
modification difference Aw,, and its effect on the weak-mixing angle measurement with 
a nuclear target. 

ACKNOWLEDGMENTS 

S.K. was supported by the Grant-in- Aid for Scientific Research from the Japanese 
Ministry of Education, Culture, Sports, Science, and Technology. 

REFERENCES 

1. G. p. Zeller et al, Phys. Rev. Lett. 88, 091802 (2002); Phys. Rev. D65, 1 1 1 103 (2002). 

2. S. Kumano, Phys. Rev. D66, 1 1 1301 (2002). 

3. M. Hirai, S. Kumano, and T.-H. Nagai, Phys. Rev. D71, 1 13007 (2005). 

4. M. Hirai, S. Kumano, and T.-H. Nagai, Phys. Rev. C70, 044905 (2004); M. Hirai, S. Kumano, and M. 
Miyama, Phys. Rev. D64, 034003 (2001). 



